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Cohesin-based chromatin interactions enable 
regulated gene expression within preexisting 
architectural compartments 
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Chromosome conformation capture approaches have shown that interphase chromatin is partitioned into spatially 
segregated Mb-sized compartments and sub-Mb-sized topological domains. This compartmentalization is thought to fa- 
cilitate the matching of genes and regulatory elements, but its precise function and mechanistic basis remain unknown. 
Cohesin controls chromosome topology to enable DNA repair and chromosome segregation in cycling cells. In addition, 
cohesin associates with active enhancers and promoters and with CTCF to form long-range interactions important for 
gene regulation. Although these findings suggest an important role for cohesin in genome organization, this role has not 
been assessed on a global scale. Unexpectedly, we find that architectural compartments are maintained in noncycling 
mouse thymocytes after genetic depletion of cohesin in vivo. Cohesin was, however, required for specific long-range 
interactions within compartments where cohesin-regulated genes reside. Cohesin depletion diminished interactions be- 
tween cohesin-bound sites, whereas alternative interactions between chromatin features associated with transcriptional 
activation and repression became more prominent, with corresponding changes in gene expression. Our findings indicate 
that cohesin-mediated long-range interactions facilitate discrete gene expression states within preexisting chromosomal 
compartments. 



[Supplemental material is available for this article.] 

The regulated transcription of mammalian genomes packaged into 
nuclei six orders of magnitude smaller than the length of chro- 
mosomal DNA requires a complex organization. The underlying 
mechanisms are beginning to be explored in terms of the bio- 
physical properties of the DNA polymer and associated chromatin 
(Fudenberg and Mirny 2012) and nuclear landmarks such as nu- 
clear lamina that provide genomic scaffolds (Guelen et al. 2008; 
McCord et al. 2013; Meuleman et al. 2013). Of particular interest 
are functional interactions that operate within these constraints to 
facilitate long-range interactions between gene regulatory ele- 
ments (Merkenschlager and Odom 2013). 

Genome-scale chromosome conformation capture has shown 
that interphase chromatin is organized into Mb-scale compart- 
ments that are "open/' gene-rich, highly transcribed, and inter- 
active (A) or "closed" (B), gene-poor, and less transcriptionally 
active (Lieberman-Aiden et al. 2009). Compartments are com- 
posed of topologically associated domains (TADs), which comprise 
hundreds of kilobases and are relatively invariant between differ- 
ent cell types (Lieberman-Aiden et al. 2009; Dixon et al. 2012; Nora 
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et al. 2012; Gibcus and Dekker 2013). This architecture is thought 
to facilitate the regulation of gene expression by constraining the 
number of regulatory elements a given gene is likely to encounter 
to those that are colocated within the same compartment or do- 
main (Gibcus and Dekker 2013). It is unknown how compartments 
and domains are built and how they contribute to the precise 
regulation of gene expression. 

The cohesin complex is essential for genome integrity in 
cycling cells, where it facilitates post-replicative DNA repair and 
sister chromatid cohesion by controlling chromosome topology 
(Nasmyth and Haering 2009; Aragon et al. 2013). In addition, 
cohesin contributes to the regulation of gene expression by 
mechanisms thought to involve long-range interactions between 
its binding sites at regulatory elements associated with CTCF 
(Parelho et al. 2008; Rubio et al. 2008; Stedman et al. 2008; Wendt 
et al. 2008) or with active promoters and enhancers (Misulovin 
et al. 2008; Kagey et al. 2010; Schmidt et al. 2010; Faure et al. 2012; 
Dorsett and Merkenschlager 2013). Taken together these proper- 
ties suggest a role for cohesin in genome organization. 
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Cohesin-based chromatin interactions 



Elucidating the global contribution of cohesin to the orga- 
nization of the genome remains a challenge, not least because the 
depletion of cohesin from proliferating cells interferes with DNA 
replication, DNA repair, and chromosome segregation (Nasmyth 
and Haering 2009). Here we use a genetic approach (Seitan et al. 
2011) to define the contribution of cohesin to the organization of 
the genome and the regulation of gene expression in noncycling 
cells in vivo. 

Contrary to expectation, we find no major role for cohesin in 
the maintenance of architectural features of genome organization 
in our experimental system. However, cohesin depletion reduced 
long-range interactions between cohesin-bound sites, and the 
resulting chromosomal interaction landscape was characterized by 
alternative interactions between chromatin features associated 
with transcriptional activation and repression. Interestingly, this 
reorganization of long-range interactions was accompanied by 
changes in gene expression in which genes at the low end of the 
expression spectrum were preferentially 
up-regulated, whereas genes at the high 
end of the expression spectrum were 
preferentially down-regulated. Our data 
indicate that the organization of the ge- 
nome into architectural compartments 
and the random assortment of genes and 
regulatory elements within them are in- 
sufficient for the precise regulation of 
gene expression. Rather, cohesin en- 
ables discrete gene expression states 
by promoting cohesin-based interac- 
tions within a preexisting architectural 
framework. 



deleted thymocytes reflected that of total RAD21 (Supplemental 
Fig. SI). ChIP and qPCR showed that RAD 21 binding to known 
cohesin sites was reduced by 80%-90% (Supplemental Fig. SI). 
This system allowed us to investigate the role of cohesin in gene 
expression and the structural and functional organization of the 
genome in noncycling cells without resorting to cell lines or in 
vitro culture systems. 

We prepared Hi-C libraries from two biological replicates of 
control and cohesin-deficient thymocytes and obtained a total 
of —392M unique valid pairs, —203M for control and —188M 
for cohesin-deficient thymocytes. To define chromosomal com- 
partments, we applied eigenvector analysis of chromosomal or- 
ganization to iteratively corrected genome-wide maps (Fig. 1A; 
Imakaev et al. 2012) at 140-kb resolution. Strikingly, the as- 
signment of chromosomal compartments was highly correlated 
for control and cohesin-deficient thymocytes (Fig. 1B,C). We 
conclude that, unexpectedly, Mb-scale architectural compart- 
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Results 

To define the contribution of cohesin 
to the organization of the genome and 
the regulation of gene expression in 
noncycling cells, we deleted the locus 
encoding the cohesin subunit RAD21 
by the activity of a CD4Cre transgene 
when developing thymocytes exit the 
cell cycle as part of their developmental 
program in vivo. Developing thymo- 
cytes arrest at the Gl phase of the cell 
cycle when DNA is unreplicated and 
chromosomes are present as single cop- 
ies, not as sister chromatids. In contrast 
to some model organisms, homologous 
chromosomes are not paired in mam- 
malian interphase, and cohesin is not 
tasked with holding sister chromatids or 
homologs together. Rad21 lox/lox CD4Cre 
CD4 + CD8 + small double positive thy- 
mocytes (referred to as "cohesin-deficient 
thymocytes" below) showed a significant 
reduction in total RAD21 protein ex- 
pression (79%-93%, depending on the 
reference used for normalization) (Sup- 
plemental Fig. SI). Most of the cohesin 
present in thymocytes was associated 
with chromatin, and the decrease in 
chromatin-associated RAD21 in Rad21- 
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Figure 1. Chromosomal compartments are resilient to the depletion of the cohesin subunit RAD21 
from noncycling thymocytes in vivo. (A) Eigenvector analysis of chromosomal organization. (B) Com- 
partment tracks (top) and interaction matrices (bottom) for chrl 3 at 140-kb resolution in control and 
cohesin-deficient thymocytes. Of 1 7,548 regions evaluated, only two showed a consistent compart- 
ment change from A to B and two from B to A (change in eigenvector value >0.03): chrl 8: 3640001 and 
chrl 8: 33600001. (C) Hi-C compartment data are highly correlated for control and cohesin-deficient 
thymocytes. 
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ments assigned by eigenvector analysis of chromosomal orga- 
nization are resilient to reduced RAD21 protein expression in 
interphase. 

Cohesin binding predicts perturbed long-range interactions 
in cohesin-deficient thymocytes 

We applied the HOMER pipeline (Heinz et al. 2010) to our Hi-C 
data (see Methods, section on Hi-C data analysis, for details). We 
first identified 100-kb genomic regions that interacted signifi- 
cantly with each other in either control or cohesin-deficient cells 
(FDR = 0.1), taking into account linear genomic distance and se- 
quencing depth. Interactions between these 100-kb regions were 
then compared between control and cohesin-deficient cells (Fig. 2A). 
Although the majority of significant interactions were unchanged 
in cohesin-deficient cells, 10,917 interactions were significantly 
altered in the pooled data, of which 1476 were found in replicate 1 
and 5004 in replicate 2 (P < 0.05). There was a highly significant 
overlap of 502 differential interactions between replicates (P < 
10~ 15 , odds ratio of 5.25 for 224 increased interactions and 9.96 for 
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278 decreased interactions) (Fig. 2B). All differential interactions 
were intra-chromosomal, and the great majority were contained 
within individual A compartments (Fig. 2C; see Supplemental 
Table 1 for a list of differential interactions determined by 
HOMER). 

To explore the relationship between cohesin binding and 
cohesin-dependent differential interactions, we focused on the 
946 distinct 100-kb genomic regions that participated in the 502 
differential interactions. These differentially interacting regions 
(Supplemental Table 1) were significantly enriched for the binding 
of the cohesin subunit RAD21, CTCF, the cohesin loading factor 
NIPBL, and the mediator subunit MED1 (Fig. 2D; Supplemental 
Fig. S1C,D for validation of the MED1 and NIPBL antibodies; ChIP 
and Supplemental Fig. SI E for an analysis of the relationship be- 
tween RAD21, CTCF, and NIPBL ChlP-seq peaks). 

At the 100-kb level, decreased interactions in particular were 
enriched for cohesin, both with and without CTCF (cohesin-non- 
CTCF or CNC in Fig. 2D), CTCF, and NIPBL (Fig. 2D). Increased 
interactions were enriched for marks of transcriptional activity, 
including MED1, H3K4me3, and RNA Pol II (Fig. 2D). The observed 
decrease in interactions between regions 
rich in cohesin and CTCF binding is in- 
tuitive based on reduced interactions be- 
tween cohesin and CTCF binding sites in 
cohesin-depleted cells (Hadjur et al. 2009; 
Kagey et al. 2010; Seitan et al. 2011). In- 
creased interactions between regions 
rich in features associated with active 
genes (MED1, H3K4me3, and RNA Pol II) 
could suggest a role for transcriptional 
activity in such interactions (Fraser and 
Bickmore 2007), although the 100-kb 
resolution achieved in the analysis of 
such regions is insufficient to pinpoint 
the role of individual features in driving 
differential interactions (see below). 

Interactions that decreased upon 
cohesin depletion were strong in con- 
trol cells, whereas interactions that in- 
creased following cohesin depletion 
were significantly weaker than average 
in control cells where cohesin was pres- 
ent (Fig. 2E). 



Cohesin depletion perturbs gene 
expression in open compartments 

RNA-seq analysis of two independent 
biological replicates at a depth of 144M 
total reads defined 1153 genes that were 
differentially expressed between control 
and cohesin-deficient thymocytes (FDR = 
0.05; 703 up-regulated, 450 down-regu- 
lated) (Fig. 3 A). We validated 15 of these 
by quantitative reverse transcriptase poly- 
merase chain reaction (qRT-PCR) and found 
that both methods were in close agree- 
ment (Fig. 3B). Comparison of RNA-seq 
data with Hi-C compartment analysis 
showed that 98% of deregulated genes 
reside in open (A) compartments (P < 
10" 15 , odds ratio = 5.49) (Fig. 3C). 
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Genes that are sensitive to cohesin dosage are bound by 
cohesin, CTCF, and NIPBL 

The great majority of genes that were deregulated in cohesin- 
deficient thymocytes were bound by RAD21, NIPBL, or CTCF (Fig. 
3D), indicating that most are direct targets of cohesin-mediated 
regulation. In agreement with this, the most highly enriched gene 
ontology (GO) terms (adjusted P< 10~ 8 ) included functions related 
to development, transcription, signal transduction, lymphocyte 
activation, and differentiation as well as hematopoiesis and the 
immune system (Fig. 3E). There was no enrichment for cell cycle 
related terms including cell cycle, DNA replication, chromosome 
segregation, and checkpoint activation. Terms relating to DNA 
damage, DNA repair, and cell division that were highly enriched in 
previous studies using RNAi-mediated depletion of cohesin from 
dividing ES cells (Kagey et al. 2010) were not enriched in cohesin- 
deficient thymocytes, validating our rationale for an experimental 
system based on nondividing cells. 

Taken together with the observation that architectural com- 
partments remain largely intact in cohesin-deficient thymocytes, 
these results suggest that deregulated gene expression occurs at the 
level of individual loci, and is not secondary to a global collapse in 
genome organization. 

Examples of deregulated genes located in differentially 
interacting regions as determined by HOMER analysis of Hi-C data 
included the region just 3' of the T cell receptor locus Tcra on chrl4 
and around the Cd3 gene cluster on chr9 (see Supplemental Table 1 
for a list of differentially expressed genes within differentially in- 
teracting regions in cohesin-deficient thymocytes). 

The Tcra locus is highly expressed in developing thymocytes, 
undergoes V(D)J recombination, and is separated from the 3' 
flanking region by CTCF and cohesin sites with CTCF-dependent 
insulator activity (Zhong and Krangel 1999; Magdinier et al. 2004). 



Figure 2. Cohesin binding predicts perturbed long-range interactions in cohesin-deficient thymo- 
cytes. (A) Schematic of approach for the mapping of differentially interacting regions. (B) Circos plot 
(Krzywinski etal. 2009) illustrating the chromosomal position of differential interactions in the context 
of chromosomal compartments. The HOMER software suite was used to determine significant in- 
teractions between 1 00-kb genomic regions in either control or cohesin-deficient thymocytes (FDR = 
0.1 ; replicates pooled). Of 1 0,91 7 interactions that were significantly altered in the pooled samples, 1 476 
interactions changed in replicate 1 and 5004 in replicate 2 (P < 0.05; the remaining interactions were only 
seen in the pooled data). Of 502 differential interactions that were shared between replicates, 278 were 
decreased (blue) and 224 were increased (red). All differential interactions were intra-chromosomal. 
Compartment assignment is indicated at the base of the interactions: black for open compartments and 
gray for closed compartments. (C) Differential interactions in cohesin-deficient cells are largely contained 
within preexisting chromosomal compartments. (A) Differential interactions entirely contained within 
the same A compartment. (B) Differential interactions entirely contained within the same B compart- 
ment (the frequency of intra-B interactions was <2% among down-regulated interactions). (A-B) Dif- 
ferential interactions bridging A and B compartments. (Other) interacting regions are either unassigned 
(i.e., inconsistent between replicates) or bridge two distinct A or B compartments. Down-regulated 
interactions (top) and up-regulated interactions (bottom) are shown separately. (D) Features enriched in 
100-kb regions that show differential interactions in cohesin-deficient thymocytes. Differential in- 
teractions between control and cohesin-deficient samples involved 946 unique 100-kb regions (510 
involved in decreased interactions, 427 in increased interactions, and 9 involved in both) that partici- 
pated in 502 unique differential interactions (278 that were decreased and 224 that were increased) and 
that were shared between replicates (P < 0.05). We tested whether differentially interacting 100-kb 
regions were enriched for the presence of Rad21 , CTCF, NIPBL, MED1 , H3K4me3, and RNAP2 binding 
events. See Supplemental Figure SI C,D for validation of the MED1 and NIPBL antibodies and ChIP, and 
Supplemental Figure SI E for an analysis of the relationship between RAD21 , CTCF, and NIPBL ChlP-seq 
peaks. Differentially interacting regions were significantly enriched for the binding of the cohesin 
subunit RAD21 , both with and without CTCF (cohesin-non-CTCF; CNC), and of the cohesin-associated 
factors CTCF and NIPBL, and features of transcriptional activity including MED1, H3K4me3, and RNAPol 
II. Differential interactions were further classified into decreased and increased interactions. (£) Strength 
distribution of cohesin-dependent interactions. Using the number of Hi-C reads as an indicator of the 
strength of interactions, differential interactions that are decreased in cohesin-depleted thymocytes are 
similar in strength to unchanged interactions in control cells, whereas increased interactions tend to be 
weak before cohesin depletion (Mann-Whitney latest P< 10~ 15 ). Outliers are not depicted. 



Depletion of cohesin results in increased interactions between Tcra 
and the 3 ' flanking region, which are visible in the Hi-C data and 
consistent with a scenario in which deregulated interactions with 
the Tcra enhancer (Ea) affect the expression of genes in the region 
(Fig. 4A). 

HOMER identified differential interactions around the highly 
expressed Cd3 gene cluster (Fig. 4B, top panel). The Bcl9 and Mpzl2 
genes are located in regions of increased interactions and are up- 
regulated in cohesin-deficient thymocytes (Fig. 4B, top panel). 
Chromosome conformation capture (3C) analysis within the dif- 
ferentially interacting 100-kb region containing the Cd3d, Cd3e, 
Cd3g, and Mpzl2 (Fig. 4B, top panel [ii]) showed that the depletion 
of cohesin significantly reduced local interactions between CTCF 
and RAD21 sites (Fig. 4B, bottom panel) for three of four inter- 
actions tested (Fig. 4B, inset, bottom left). Reduced local inter- 
actions may facilitate the increased long-range interactions detected 
at the 100-kb level and could also be relevant to altered gene ex- 
pression in this region: CTCF sites demarcate the Cd3 gene cluster 
from the neighboring Mpzl2 gene that is barely expressed in con- 
trol thymocytes, but up-regulated in cohesin-deficient thymocytes 
(RNA-seq data in Figure 4B, bottom panel; and RT-PCR validation, 
inset, bottom right). 

Our data causally link cohesin to long-range chromatin in- 
teractions. Although the associated changes in gene expression 
may be rationalized by reduced boundary or insulator function, 
the link between long-range interactions and deregulated gene 
expression remains correlative at this time. 

Predictive features of genes that show cohesin-dependent 
expression 

To delineate factors associated with gene expression changes in 
cohesin-deficient thymocytes, we used a regression model that 
integrated gene expression, Hi-C, and 

ChlP-seq data. We assigned genes to one 

of three classes (up-regulated, down-reg- 
ulated, or unchanged) and tested for the 
presence or absence of ChlP-seq peaks 
near each gene promoter (TSS ±2.5 kb) as 
well as gene location within 100-kb re- 
gions that interact differentially in con- 
trol and cohesin-deficient thymocytes. 
Differential interactions were further 
divided into interactions that were stron- 
ger (DI region, Up) or weaker (DI region, 
Down) in cohesin-deficient thymocytes. 
We also considered the presence of the 
H3K4me3 histone modification; the bind- 
ing of RAD21, CTCF, NIPBL, and Mediator, 
RNA Pol II, paused RNA Pol II at the pro- 
moter (Hendrix et al. 2008); the presence 
of promoter CpG islands (CGI); and gene 
length. We determined the relative im- 
portance of each variable in the multi- 
variate model and ranked the corre- 
sponding coefficients by their statistical 
significance. The presence of a CpG 
island (CGI) at the promoter, gene length, 
the presence near promoters of RAD21, 
particularly without CTCF (CNC), and 
promoter-associated RNA Pol II emerged 
as important variables (Fig. 4C; see Sup- 
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Figure 3. Gene expression in cohesin-deficient thymocytes. (A) Volcano plot of RNA-seq data. Of the 1 7,849 genes assayed, 1153 are significantly 
differentially expressed in cohesin-deficient thymocytes (FDR = 0.05), 450 of which are up-regulated (orange) and 703 down-regulated (green). The red 
dot represents Rad21 . Eight genes above the -logl O(P-value) = 50 threshold (all up-regulated) were omitted from the plot. (B) Validation of RNA-seq data 
for 15 transcripts over a wide range of expression levels. Shown is fold-change measured by qRT-PCR (x-axis, two independent biological replicates) 
against the fold-change measured by RNA-seq (/-axis) and the Pearson's correlation coefficient (r) for two independent RNA-seq experiments. (C) Genes 
affected by cohesin depletion are associated with open compartments. Approximately 98% of differentially expressed genes that could be assigned to 
compartments (1121 of 1142 autosomal and X-linked genes) reside in open (A) compartments, compared with 91.2% (16,255 of 17,819) of genes 
included in our analysis. Genes spanning more than one compartment were assigned to A when they overlapped at least partially with open compart- 
ments; genes overlapping compartments that could not be clearly defined as A or B because they were inconsistent between replicates were called 
unassigned. Only 21 deregulated genes were found outside open compartments. Of these, 1 1 were located in unassigned compartments and just 1 0 of 
1142 deregulated genes were located in closed (B) compartments. This corresponds to a highly significant depletion of deregulated genes in B com- 
partments (P < 1 0 , odds ratio = 0.1 2). (D) Cohesin-regulated genes are bound by cohesin and associated factors. Associations are shown for ChlP-seq 
peaks for RAD21 , NIPBL, and CTCF (Shih et al. 201 2) within 2.5 kb of transcription start sites (TSS) and within 1 0 kb of canonical gene bodies. (£) Gene 
Ontology (GO) analysis of genes that are differentially expressed in cohesin-deficient thymocytes. Representative GO biological process terms with 
adjusted P< 10~ 8 are shown (the complete list can be found in Supplemental Table 2). 



plemental Fig. S2 for the results of a univariate analysis considering 
each variable separately). Interestingly, location within regions 
that showed differential interactions in cohesin-deficient thymo- 
cytes — but not a control set of interacting regions that showed 
no differences between control and cohesin-deficient thymocytes 
("Random DI region" in Fig. 4C) — was predictive of gene expression 
changes. Decreased gene expression in particular was associated with 
differential interactions (P-value < 10~ n , odds ratio = 2.63) (Fig. 4C). 

Investigating variables that were highlighted by the regres- 
sion model, we first probed the relationship between deregulated 
gene expression and differential long-range interactions. This 
analysis showed a strong association between decreased interac- 
tions and down-regulated gene expression (Supplemental Fig. S3 A). 
Further analysis confirmed that genes with increased promoter CpG 
density were more likely to be up- than down-regulated in cohesin- 
deficient thymocytes (Supplemental Fig. S3B, left), and also more 
likely to coincide with differentially interacting regions in cohesin- 
deficient cells (Supplemental Fig. S3B, right). Moreover, longer 



genes were more likely to be down-regulated in cohesin-deficient 
cells (Supplemental Fig. S3C, left) and were found preferentially 
in regions of reduced interactions (Supplemental Fig. S3C, right). 
A set of complex loci of significantly greater length than average 
(Mann-Whitney [/-test P < 10~ 15 ) overlap ultraconserved non- 
coding elements, regulatory blocks or archipelagos (Akalin et al. 
2009; Montavon and Duboule 2012; Dimitrieva and Bucher 2013). 
These loci were preferentially down-regulated in cohesin-deficient 
thymocytes (Supplemental Fig. S3D), which is consistent with 
a role for cohesin in long-range regulation of genes with complex 
regulatory inputs. 

Cohesin depletion perturbs long-range interactions 
within architectural compartments and compresses 
the dynamic range of gene expression 

We used Structured Interaction Matrix Analysis (SIMA) (Lin et al. 
2012) to obtain a high-resolution view of interactions between 
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Figure 4. Features and predictors of cohesin-dependent gene expression. (A) The T cell receptor locus Tcra and 3' flanking region of chrl 4 are shown 
alongside Hi-C heat maps at 1 40-kb resolution for control (top) and cohesin-deficient thymocytes (bottom), epigenomic features, and genes that are up- 
(red) or down-regulated (blue, no examples in this region). Differentially interacting (Dl) regions as determined by HOMER analysis of Hi-C data are shown 
in red (up-regulated) or blue (down-regulated). Dashed gray lines demarcate increased interactions with the region 3' of Tcra. (B) HOMER-identified 
increased and decreased interactions around the Cd3 gene cluster on chr9. (Top) increased interactions marked as (i) chr9:44300000-44400000 and 
(ii) chr9:44800000-44900000 contain the Bcl9 and Mpzl2 genes that are up-regulated in cohesin-deficient thymocytes. (Bottom) Gene expression and 3C 
analysis of region (ii). Hind III fragments containing CTCF and RAD21 binding sites are shaded gray. (Left inset) Primary 3C data (mean ± SD, n = 3). (Right 
inset) RT-PCR validation of Mpzl2 expression (mean ± SD, n= 3). (C) Multinomial logistic regression model integrating gene expression, Hi-C, and ChlP- 
seq data to predict up-regulated, down-regulated, and unchanged genes. We tested each gene for the following features: (1 ) the presence or absence of 
ChlP-seq peaks near the promoter (TSS ± 2.5 kb); (2) location within 1 00-kb differentially interacting regions (divided into interactions that were stronger, 
Dl region [Up], or weaker, Dl region [Down], in cohesin-deficient thymocytes); (3) the presence of the H3K4me3 histone modification; (4) the binding of 
RAD21 , CTCF, NIPBL and Mediator (MED1 ), RNA Pol II (RNAP2), paused RNA Pol II at the promoter (Hendrix et al. 2008); (5) the presence of a promoter 
CpG island (CGI); and (6) gene length. Error bars represent 95% confidence intervals. Variables are ranked by coefficient significance from left to right. Of 
17,849 genes assayed by RNA-seq, 1461 resided in Dl regions. Of 450 down-regulated genes, 83 resided in Dl regions, which represents a strong 
enrichment (P-value < 1 0 -11 , odds ratio 2.63). Up-regulated genes were only slightly enriched in Dl regions (P < 0.05, odds ratio 1 .27). 
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specific chromatin features within the chromosomal compart- 
ments assigned by eigenvector analysis of chromosomal organi- 
zation (Fig. 5 A). Although the number of Hi-C reads that map to 
any one specific feature — for example a cohesin-bound site — is too 
low to assign interactions with confidence, this approach is de- 
signed to combine information for multiple occurrences of this 
feature (Lin et al. 2012). We focused on open (A) compartments, 



where the great majority of cohesin binding sites and cohesin- 
regulated genes reside (Fig. 3C), and selected a range of features for 
analysis. These included ChlP-seq peaks for cohesin and associated 
factors, histone modifications indicative of active (H3K4me3), and 
repressed (H3K27me3) chromatin states, transcription start sites of 
active and silent genes, as well as control sites that did not overlap 
with these features. For each compartment (1, 2,. . ., n in Fig. 5 A) we 
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Figure 5. (Legend on next page) 
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counted Hi-C reads connecting features of the same type (homo- 
typic interactions, illustrated for feature 1 in Fig. 5A) and Hi-C 
reads connecting different features (heterotypic interactions, il- 
lustrated for features 2 and 3 in Fig. 5 A), associating Hi-C reads that 
mapped within 10 kb of each feature (Lin et al. 2012). To determine 
the impact of cohesin on these interactions we compared inter- 
actions in control and cohesin-deficient cells for each feature (and 
pair of features) within open compartments (Methods; Supple- 
mental Table 5). 

This approach revealed that interactions between RAD21 and 
CTCF sites were reduced in cohesin-deficient thymocytes (Fig. 5B), 
and the analysis of interactions between all pairs of features 
revealed that RAD21-RAD21 and RAD21-CTCF interactions were 
most strongly decreased (Fig. 5C). This result is consistent with 
previous 3C experiments that indicated reduced interactions be- 
tween such sites in cohesin-depleted cells (Fig. 4B; Hadjur et al. 
2009; Kagey et al. 2010; Seitan et al. 2011). 

Interactions that remained — or even increased — in cohesin- 
deficient cells included features of active transcription such as 
NIPBL, H3K4me3, and RNA Pol II and extended to interactions 
between activation-associated features and repressive H3K27me3 
marks (Fig. 5B,C). Upon cohesin depletion, promoters of genes 
that were silent in control cells also showed detectable interac- 
tions with a range of features linked to transcriptional activation 
(Fig. 5C). 

To determine whether these interactions in cohesin-deficient 
cells were selective, we identified control sites that were at least 
10 kb removed from other features. Interactions involving these 
"random" sites showed little cohesin dependence, suggesting that 
increased interactions preferentially involve sites that are marked 
by the features analyzed, and thereby provide a measure of selec- 
tivity (Supplemental Fig. S4). 

Since SIMA analysis is based on averaged occurrences of 
the specific features we do not know whether the increased rep- 
resentation of alternative interactions in cohesin-deficient cells is 
absolute or relative to the loss of cohesin-based interactions. Ei- 
ther way, the data indicate a shift in the chromatin landscape in 
cohesin-deficient cells from cohesin-based to alternative interac- 
tions (Fig. 5D). 

As a proxy for the scale of differential interactions with de- 
creased and increased representation in our Hi-C data, we stratified 
SIMA results by compartment size. Interestingly, the reduction 



of cohesin-based interactions was most pronounced in compart- 
ments <1 Mb, suggesting an upper limit for cohesin-based inter- 
actions within compartments (Fig. 5E; Supplemental Fig. S5). The 
alternative interactions detected in cohesin-deficient cells in- 
creased with compartment sizes >1 Mb, suggesting that they 
preferentially occur over larger distances or in larger compart- 
ments, which may have a higher degree of complexity (Fig. 5E; 
Supplemental Fig. S5). 

The prominence of interactions between features associated 
with both active gene expression and silencing in cohesin-de- 
ficient cells (Fig. 5C) resonates with models where cohesin not 
only facilitates specific interactions, but also provides separation 
between genes and regulatory elements (Wallace and Felsenfeld 
2007; Phillips and Corces 2009; van Steensel 2011; Dixon et al. 
2012). To address how the altered chromatin interactions in cohesin- 
deficient cells may impact on gene expression, we stratified genes 
according to their level of expression. Gene expression was per- 
turbed across the entire range of the expression spectrum (Fig. 6A). 
In addition, genes with low expression were more often up- 
regulated, whereas genes with high expression were more often 
down-regulated (Fig. 6A). This resulted in a systematic skewing of 
gene expression away from the extremes of the dynamic range and 
toward average values (Pearson's correlation coefficient r= -0.16, 
P < 10~ 7 ) (Fig. 6A). Accordingly, the proportion of up-regulated 
genes was higher at the low end of the expression spectrum, 
whereas the proportion of down-regulated genes was higher at the 
high end of the expression spectrum (Pearson's correlation co- 
efficient r- -0.97, P < 10~ 4 ) (Fig. 6B). Hence, genes showed a more 
uniform expression in cohesin-deficient cells. 

Discussion 

The cohesin complex provides physical linkage between sister 
chromatids from the time of chromosome duplication in S-phase 
until chromosomes segregate in cell division (Nasmyth and Haering 
2009) and can form long-range interactions that link gene regu- 
latory elements with their targets in interphase (Hadjur et al. 2009; 
Kagey et al. 2010; Seitan et al. 2011). Unexpectedly, our Hi-C 
analysis shows that the maintenance of compartments in non- 
dividing mammalian cells was not affected by the depletion of 
cohesin, at least not at a level we can achieve in our experimental 
system. This reduction is far in excess of the 20%-30% reduction 



Figure 5. The chromatin landscape of cohesin-deficient thymocytes is characterized by the loss of cohesin-based interactions and the detection of 
alternative interactions. (A) Outline of approach used to score interactions between specific chromatin features located within architectural compartments 
as assigned by eigenvector analysis of chromosomal organization (SIMA) (Lin et al. 201 2). Open compartments (1 , 2, . . ., n) are indicated by gray boxes. 
Orange, green, and black symbols (features 1 , 2, etc.) represent chromatin features, and purple lines indicate "featureless" control sites. Hi-C interactions 
between features are indicated by dashed lines. For each compartment, we counted Hi-C reads connecting features of the same type (homotypic 
interactions, illustrated for feature 1), and Hi-C reads connecting different features (heterotypic interactions, illustrated for features 2 and 3), assigning Hi-C 
reads that mapped within 1 0 kb of each feature (Lin et al. 201 2). We compared interactions in control and cohesin-deficient cells for each feature and pair 
of features within each open compartment and combined these scores as a measure for the cohesin dependence of long-range interactions between 
specific features. (B) Impact of cohesin deficiency on homotypic cohesin-based and alternative interactions. SIMA was used to determine the enrichment of 
Hi-C reads in interactions connecting "like" features in control and cohesin-deficient thymocytes. The difference between normalized enrichment ratios in 
each condition was assessed by the Wilcoxon signed-rank test (see Methods). (C) Impact of cohesin deficiency on all pairwise feature-based interactions. 
SIMA results for homotypic interactions are shown together with those associated with interactions connecting different features (heterotypic in- 
teractions). Refer to B for details. Homotypic RAD21-RAD21 interactions are the most decreased, followed by CTCF-RAD21 interactions. Interactions 
between features associated with active transcription are strongly increased in cohesin-deficient thymocytes, as are interactions involving the repressive 
histone modification H3K27me3 (Zhang et al. 201 2) with marks of active transcription. (D) Cytoscape representation of SIMA results in C Edge color and 
width correspond to the Wilcoxon signed-rank test effect size and significance, respectively. (£) Length scale of lost and gained cohesin-dependent feature- 
based interactions. Boxplots are based on SIMA comparisons stratified into three classes according to compartment size (<1 Mb, 1 -3 Mb, >3 Mb). Effect 
sizes for decreased interactions (RAD21 -RAD21, CTCF-RAD21, CTCF-CTCF) and increased interactions (remainder) were grouped and are indicated 
separately. The effect of reduced cohesin-based interactions is most pronounced within smaller compartments and decreases when larger compartments 
are considered (Pearson's correlation coefficient r= -0.73, P< 0.05), whereas alternative interactions increased with compartment size (r= 0.7, P< 0.001 ). 
Outliers are not depicted. 
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Figure 6. Cohesin depletion compresses the dynamic range of gene 
expression. (A) Genes were stratified into 10 equally sized log intervals 
from low (0-1) to high (>9) based on the average gene expression of 
control and cohesin-deficient thymocytes. Boxplots indicate the distri- 
bution of gene expression fold changes in cohesin-deficient thymocytes. 
The number of genes in each bin is indicated (bins 1 and 2 are empty). 
Note that lowly expressed genes are frequently up-regulated, whereas 
highly transcribed genes tend to be down-regulated. Note that this pat- 
tern does not result from ascertainment bias in which lowly expressed 
genes can only be up-regulated and vice versa, because we stratified 
genes according to the mean of their expression in control and cohesin- 
depleted cells. Therefore, the direction of regulation is the inverse when 
control cells are compared to cohesin-deficient cells (not shown). P-values 
are based on one-sample Wilcoxon signed-rank tests and indicate signif- 
icant difference from zero (no change). Genes with zero mean expression 
in both cohesin-deficient and control thymocytes are excluded and out- 
liers are not depicted. (B) Bar plot indicating the proportion of up- and 
down-regulated genes in each gene expression interval (see A). The pro- 
portion of up-regulated genes is anti-correlated with interval rank 
(Pearson's correlation coefficient r= -0.97, P< 10 4 ). 



known to cause severe developmental abnormalities in model or- 
ganisms and in human Cornelia de Lange syndrome (CdLS), sug- 
gesting that a breakdown of genome organization at the level of 
compartments may not be the cause for deregulated gene expres- 
sion in CdLS (Strachan 2005; Kawauchi et al. 2009). In contrast to 
the preservation of architectural genome organization, we find 
that cohesin depletion alters long-range chromosomal interac- 



tions within compartments and results in a more uniform ex- 
pression of genes affected by cohesin depletion. 

We found that the density of cohesin binding sites correlated 
with differential interactions of 100-kb regions that showed both 
increased and decreased interactions in cohesin-deficient thymo- 
cytes. Analysis at the local level showed that cohesin depletion 
reduced Hi-C reads that mapped within 10 kb of RAD21 and CTCF 
binding sites. Other interactions remained or even increased in 
cohesin-deficient cells. These alternative interactions included 
features of transcriptionally active and repressed chromatin. Using 
compartment size as a proxy for the length scale of interactions, 
cohesin-based interactions and the alternative interactions de- 
tected in cohesin-deficient cells appeared to differ in scale. Con- 
sistent with correlative data (de Wit et al. 2013; Phillips-Cremins 
et al. 2013), cohesin-dependent interactions within compartments 
were mostly confined to <1 Mb, i.e., the scale of topologically as- 
sociated domains (TADs). In contrast, the alternative interactions 
detected in cohesin-deficient cells increased with compartment 
sizes >1 Mb. This result suggests that TAD-scale cohesin-based 
interactions prevent the detection of longer-range alternative 
interactions. Hence, disrupting cohesin-based interactions may 
lead to a degree of mixing between chromosomal domains beyond 
the TAD scale (>1 Mb). 

Since SIMA analysis is based on averaged occurrences of the 
specific features, we do not know whether the increases are abso- 
lute or relative to the loss of cohesin-based interactions. Either way, 
the data indicate a shift in the chromatin landscape in cohesin- 
deficient cells from cohesin-based to alternative interactions. Im- 
portantly, the homogenization of interactions that results from 
this shift is consistent with the pattern of deregulated gene ex- 
pression in cohesin-deficient thymocytes. Similar to previous ob- 
servations in Drosophila (Schaaf et al. 2009), genes at the low end of 
the expression spectrum were preferentially up-regulated, whereas 
genes at the high end of the expression spectrum were preferen- 
tially down-regulated. Perhaps this is the genome-scale equivalent 
of the ability of CTCF and CTCF-associated cohesin to mediate 
insulator and boundary functions. Although such functions were 
predicted based on correlative analysis of histone modifications 
(Cuddapah et al. 2009), lamin-associated domains (Guelen et al. 
2008; McCord et al. 2013; Meuleman et al. 2013) and long-range 
interactions maps (Dixon et al. 2012), their impact on gene ex- 
pression had been previously documented only in reporter assays 
(Wallace and Felsenfeld 2007; Phillips and Corces 2009) and at 
individual loci, based on the manipulation of individual CTCF 
sites (Guo et al. 2011) or the deletion of CTCF (Ribeiro de Almeida 
et al. 2011). We suggest that cohesin-based interactions limit the 
extent of alternative interactions to enable discrete gene expres- 
sion states. 

Our conclusion that cohesin contributes to functional in- 
teractions within preexisting chromosomal compartments con- 
trasts with cohesin's structural role in providing stable cohesion 
between sister chromatids (Nasmyth and Haering 2009). On closer 
inspection, however, this conclusion is consistent with a sub- 
stantial body of existing knowledge on how cohesin interactions 
with chromatin are regulated during the cell cycle. In the Gl phase 
of the cell cycle, chromatin-associated cohesin shows relatively 
rapid turnover and has a half-life of minutes (Gerlich et al. 2006; 
McNairn and Gerton 2009; Gause et al. 2010; Onn and Koshland 
201 1). It is only during S-phase that a subset of cohesin complexes 
is stabilized by the acetylation of specific lysine residues in the 
SMC3 subunit, which displace the cohesin unloading factor WAPL 
and allow association of the cohesin-stabilizing factor sororin to 
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extend the half-life of a subset of cohesin complexes (Gerlich et al. 
2006; Kueng et al. 2006; Gause et al. 2010). This subset of long- 
lived cohesin complexes is thought to have a structural role by 
mediating cohesin between sister chromatids from S-phase until 
the transition between metaphase and anaphase, which can be 
hours or even decades in the case of human oocytes (Jessberger 
2012). However, a population of cohesin complexes with high 
turnover remains even after DNA replication, and perhaps it is 
those complexes that continue to contribute to the regulation of 
gene expression during the S- and G2-phases of the cell cycle. The 
depletion of the cohesin removal factor WAPL demonstrates the 
consequences of rendering all cohesin complexes stable (Kueng 
et al. 2006; Tedeschi et al. 2013). Interestingly, loss of WAPL causes 
dramatic changes in the structure of interphase chromatin and 
perturbs the regulation of gene expression (Kueng et al. 2006; 
Tedeschi et al. 2013). It therefore appears that dynamic cohesin 
turnover is essential for regulated gene expression, and we specu- 
late that this behavior is reflected in the contribution of cohesin to 
genome organization in interphase that is described by our data. 

Methods 

Experimental procedures 

The conditional Rad21 allele crossed to CD4Cre and methods for 
RT- and genomic PCR, chromosome conformation capture and 
ChlP-seq have been described (Seitan et al. 2011). ChIP was per- 
formed using Abeam ab992 rabbit polyclonal antibody to RAD21 
(Seitan et al. 2011), Bethyl Laboratories A300-793A rabbit poly- 
clonal antibody to MED1 (Kagey et al. 2010), and Bethyl Labora- 
tories A301-779A rabbit polyclonal antibody to NIPBL (Kagey et al. 
2010). ChIP primers (5-3') and amplicon positions (mm9) were: 
X:100692 (ChrX:100691761-100691894): fwd: TGGCTTGTACT 
TCCAGATCAT, rev: AATGTATAAGCAAGCTAGTACGCA; 17:13025 
(mm9 Chrl7:13025211-13025311): fwd: GGTCTCGCCAGCTT 
GCTATTTC, rev: TGGCTAGATGTCATTGTGGTGG; 17:17783 
(Chrl7:17782844-17782944): fwd: AGGACATAGTCGCTTGAGTGAT 
GG, rev: TTGGGCCAGGCTGGTACTTT; 11:53442 (Chrll:53441557- 
53441657): fwd: CAGGATGAGGATATGACTAGCTGTGG, rev: TC 
CAGTCCAGTTCCAAAAGGAGC. RNA-seq was done as described 
(Seitan et al. 2011) except that we used TruSeq kits according to 
the manufacturers' instructions (Illumina). Hi-C libraries were 
prepared as described (Lieberman-Aiden et al. 2009; Belton et al. 
2012). 

Computational methods 

ChlP-seq read mapping and peak calling 

Raw read alignment, filtering, and peak-calling for RAD21, MED1, 
NIPBL CTCF, and definition of CNCs was done as previously de- 
scribed (Faure et al. 2012). ChlP-seq data for H3K27me3 (Zhang 
et al. 2012) was similarly processed, except CCAT version 3.0 (Xu 
et al. 2010) was used to identify the relatively broad regions oc- 
cupied by this mark (precompiled histone modification con- 
figuration). Single base pair summit positions for H3K4me3 and 
RNAP2 obtained from the Mouse ENCODE Project (Shen et al. 
2012) were extended to a width of 200 bp to define peak regions. 

RNA-seq data analysis 

Raw reads for each condition and replicate were independently 
aligned to mouse transcript sequences (cDNA sequences from 
Ensembl version 66, NCBI37/mm9) (Flicek et al. 2013) using 
Bowtie version 0.12.8 (Langmead et al. 2009) with default pa- 



rameters. Gene expression estimates and normalized count equiv- 
alents were obtained using MMSEQ version 0.11.2 (Tuno et al. 

2011) . We used the Bioconductor R package DESeq version 1.6.1 
(Anders and Huber 2010) to determine significantly differen- 
tially expressed genes in cohesin-deficient thymocytes versus con- 
trol cells (FDR = 0.05). Empirical gene expression dispersion values 
were estimated in a condition-specific fashion (method = "per- 
condition") and used to fit a dispersion-mean relationship, where 
only the fitted values were used (sharingMode = "fit-only"). T cell 
receptor gene segments, pseudogenes, ribosomal genes, and genes 
with an aggregate exon mapability score in the lowest ten per- 
centile (UCSC Genome Browser track "wgEncodeCrgMapability- 
Align50mer") were excluded from the analysis. Expressed genes 
were defined as those having log(expression_level + 1) > 1 in con- 
trol cells; otherwise genes were considered silent. 

Hi-C data analysis 

Iterative error correction of Hi-C data was performed as described 
(Imakaev et al. 2012). The HOMER Hi-C software analysis pipeline 
(http://biowhat.ucsd.edu/homer/interactions/) was used to de- 
termine significant interactions, differential interactions and to 
perform Structured Interaction Matrix Analysis (SIMA) (Lin et al. 

2012) . Briefly, paired-end reads were trimmed to remove sequence 
following the canonical Hindlll ligation junction sequence (lbp 
mismatch allowed to account for potential star activity). Trimmed 
reads were aligned independently to the mouse reference genome 
assembly (NCBI37/mm9) using BWA (Li and Durbin 2009). Paired- 
end reads were merged and filtered to remove duplicate read pairs 
("-tbp 1"), paired-end reads likely representing continuous geno- 
mic fragments or religation events ("-removePEbg"), self-ligations 
("-removeSelfLigation"), and reads originating from regions with 
unusually high tag density ("-removeSpikes 10000 5"). Addition- 
ally, only read-pairs where both ends mapped near restriction sites 
were retained ("-both"). 

To identify differential interactions, we first determined a 
"universe" of interactions on which to focus the analysis, defined 
as significant interactions between 100-kb genomic regions in ei- 
ther control or cohesin-deficient thymocytes (replicates pooled; 
FDR = 0.01). HOMER uses the binomial distribution to determine 
significant deviations above the expected number of Hi-C reads 
occurring between two loci, where the background model takes 
into account linear genomic distance and sequencing depth. For 
each replicate, high scoring differential interactions (P< 0.05) were 
then determined by comparing Hi-C read levels between control 
and cohesin-deficient thymocytes within this subset of all possible 
interactions. Only differential interactions consistently identified 
in both Hi-C replicates were retained for downstream analysis. 

To determine genomic features associated with chromatin 
interactions, we used a method that pools Hi-C information as- 
sociated with a given set of genomic regions within a specified set 
of domains (SIMA) (Lin et al. 2012). We used default resolution 
("-res 2500") and optimal Hi-C interaction search space parameters 
("-superRes 10000"). Domains of interest were defined as merged 
adjacent 140-kb regions within open compartments at least 500 kb 
("-minDsize 500000") and not more than 5 Mb in length. Within- 
domain associations were assessed independently in control and 
cohesin-deficient thymocytes for all peak sets (RAD21, MED1, 
NIPBL, CTCF, H3K4me3, H3K27me3, RNAP2) as well as all ca- 
nonical TSSs (excluding pseudogenes; Ensembl version 66; Flicek 
et al. 2013), promoters of silent genes, expressed genes, and sig- 
nificantly differentially expressed genes. Normalizing by the 
number of expected Hi-C reads under the background model and 
comparing to the randomized average (after shuffling feature po- 
sitions 10,000 times), we obtained an enrichment ratio for each 
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genomic region, indicating the association of each feature (or 
feature pair) with interactions contained within that region. To 
determine the impact of cohesin on these interaction associations, 
we compared observed/randomized enrichment ratios in control 
and cohesin-deficient cells for each feature (or feature pair) within 
each compartment. The size and direction of change in these ratios 
in cohesin-depleted thymocytes were compared using a paired 
statistical test (Wilcoxon signed-rank test) to provide a measure for 
the cohesin dependence of long-range interactions between spe- 
cific features. See Supplemental Figure S5 for SIMA results stratified 
by compartment size. 

Multinomial logistic regression model 

We used gene features derived from ChlP-seq, RNA-seq, Hi-C, ge- 
nomic sequence, and annotation to predict gene expression 
changes in cohesin-deficient thymocytes using a multinomial lo- 
gistic regression model (nnet R package version 7.3-1) (Venables 
and Ripley 2002). Binary variables used included RAD21 (over- 
lapping CTCF), cohesin-non-CTCF (CNC), CTCF, NIPBL, MED1, 
H3K4me3 and RNAP2 peak, and CpG island (CGI) (Illingworth et al. 
2010) overlap within the gene promoter (TSS ± 2.5kb), paused 
promoter RNAP2 (stalling index >4) (Hendrix et al. 2008), gene 
body overlap with 100-kb differentially interacting (DI) regions (Up, 
increased; Down, decreased) as well as a randomly selected "control" 
set of interacting regions. We also included total gene length as 
a continuous variable. The three-class categorical response variable 
was encoded as follows: "0" nondifferentially expressed, "— 1" sig- 
nificantly down-regulated, "I" significantly up-regulated. 

Data access 

ChlP-seq, RNA-seq, and Hi-C data from this study have been 
submitted to the NCBI Gene Expression Omnibus (GEO; http:// 
www.ncbi.nlm.nih.gov/geo/) under accession number GSE48763. 
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